Long time correlations in Lagrangian dynamics: a key to intermittency in turbulence 
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New aspects of turbulence are uncovered if one considers flow motion from the perspective of a fluid 
particle (known as the Lagrangian approach) rather than in terms of a velocity field (the Eulerian 
viewpoint). Using a new experimental technique, based on the scattering of ultrasounds, we have 
obtained a direct measurement of particle velocities, resolved at all scales, in a fully turbulent flow. It 
enables us to approach intermittency in turbulence from a dynamical point of view and to analyze the 
Lagrangian velocity fluctuations in the framework of random walks. We find experimentally that the 
elementary steps in the walk have random uncorrelated directions but a magnitude that is extremely 
long-range correlated in time. Theoretically, we study a Langevin equation that incorporates these 
features and we show that the resulting dynamics accounts for the observed one- and two-point 
statistical properties of the Lagrangian velocity fluctuations. Our approach connects the intermittent 
statistical nature of turbulence to the dynamics of the flow. 
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Traditional experimental studies of velocity fluctua- 
tions in turbulence rely on velocimetry measurement at 
a fixed point in space. A local velocity probe yields time 
traces of the velocity fluctuations which are then related 
to spatial velocity profiles using the Taylor hypothesis [Q . 
In this case, the flow is analyzed in terms of the Eu- 
lerian velocity field u(x, t). One of the most peculiar 
feature of homogeneous three-dimensional turbulence is 
its intermittency, well established in the Eulerian frame- 
work H . The statistical properties of the flow depend on 
the length scale at which it is analyzed. For instance, the 
functional form of the probability of measuring an Eule- 
rian velocity increment A s u(x) — u[x + s) — u(x) varies 
with the magnitude of the length scale s. Many studies 
devoted to the understanding of this feature have been 
developed along the lines of Kolmogorov and Obhukov 
1962 pioneering ideas j3|. In this case, intermittency is 
analyzed in terms of the anomalous scaling of the mo- 
ments of the velocity increments in space. It is attributed 
to the inhomogeneity in space of the turbulent activity 
and often analyzed in terms of ad-hoc multiplicative cas- 
cade models g. Although very successful at describing 
the data, these models have failed to connect intermit- 
tency with the dynamical equations that govern the mo- 
tion of the fluid. Here, we adopt a Lagrangian point of 
view. It is a natural framework for mixing and trans- 
port problems in turbulence In addition it has been 
shown in the passive scalar problem that intermittency 
is strongly connected to the particular properties of La- 
grangian trajectories [js], In the Lagrangian approach, 
the flow is parameterized by v(%a,t), the velocity of a 
fluid particle initially at position xq. Experimentally, we 
follow the motion of a single tracer particle and we con- 
sider the increments in time of its velocity fluctuations: 
A T v(t) — v(t+r)—v(t). Our first observations |7| have es- 
tablished and described intermittency in this Lagrangian 
framework. In order to understand its origin, we pro- 



pose here a new point of view. Since our measurements 
give access to the individual motion of fluid particles, 
we study intermittency from a dynamical point of view. 
We show that the anomalous scaling in the Lagrangian 
velocity increments traces back to the existence of long- 
time correlations in the particle accelerations, i.e. the 
hydrodynamic forces that drive the particle motion. 

In order to study the motion of Lagrangian tracers, 
we need to resolve their velocity fluctuations across a 
wide range of scales. To this end, we use a confined flow 
with no mean advection, so that fluid particles remain 
for long times inside a given measurement volume. The 
tracking of small tracer particles is achieved using a new 
acoustic technique based on the principle of a "continu- 
ous Doppler sonar" . The flow volume is continuously in- 
sonified with a monochromatic ultrasound which is then 
scattered by the tracer particle j?J . This scattered sound 
is detected by two transducer arrays which yield a mea- 
surement of both the particle position, by direct triangu- 
lation, and of its velocity, from the Doppler shift. Indeed, 
for an incoming sound with frequency /o , the scattered 
sound at the receiver has frequency f(t) = fo + k.v(t) 
, where v(i) is the velocity of the tracer particle and k 
is the scattering wave vector. This frequency modula- 
tion in the acoustic signal is extracted numerically, us- 
ing a high-resolution parametric method || . Figure 1 (a) 
shows the experimental set-up and an example of a par- 
ticle trajectory; Figure 1(b) gives an example of the time 
variation of one component of its velocity. A water flow 
of the von Karman swirling type ^| is generated in- 
side a cylinder by counter-rotation at Q = 7 Hz of two 
discs with radius R = 9.5 cm, fitted with eight blades 
of height 0.5 cm and set 18 cm apart. The flow power 
consumption is e = 25 W/kg, with velocity fluctuations 
Urms = 0.98 m/s. The characteristic size of the velocity 
gradients is £ = (15u 2 m;5 / e) 1 / 2 = 880 fim, larger than 
the diameter (250 /im) of the neutrally buoyant tracer 
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particle (density of 1.06). The turbulent Reynolds num- 
ber of the flow is Re = u rms i/y = 740. The large scale 
flow is axisymmetric and the fluctuations in its center ap- 
proximate well the conditions of local homogeneous and 
isotropic turbulence. The flow is insonified at 2.5 MHz, 
with the transducers located at the flow wall. The re- 
ceiver arrays are placed at 45 degrees on each side of the 
emission direction. The measurement region is the inter- 
section of the emission and detection cones. The particles 
act as Lagrangian tracers for times longer than 1 ms (be- 
low which inertia cuts off their response) , up to times as 
long as they will stay confined insided the measurement 
volume, i.e. between one and ten Tl, the Lagrangian in- 
tegral time scale (computed from the Lagrangian velocity 
autocorrelation function) . 4000 such events are analyzed, 
for a total of 1.9 x 10 6 data points sampled at 6.5 kHz. 
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FIG. 1: (a) Experimental set-up of the von Karman flow, with 
an example of a particle 3D trajectory, 206 ms (= 9.2 Tl); 
(b) corresponding velocity variation (one component shown). 

The probability density functions (PDFs) of the La- 
grangian time increments n r (Av), are shown in Fig. 2. 
They are Gaussian at integral scale (r > Tl) and vary 
continuously towards the development of stretched ex- 
ponential tails as the time increments decrease towards 
the dissipative time scale (?], [Io|. The outer curve in 
Fig. 2 is the PDFs of Lagrangian acceleration measured 
by LaPorta et al. ^ in the same flow geometry and at 
a comparable turbulent Reynolds number. This evolu- 
tion of the PDFs n r (Ai>) leads to an anomalous scal- 
ing of the velocity structure functions (|A T t!| 9 ) ~ r^ 9 ', 
with £(q) a non-linear function of q. The intermit- 
tency spectrum is usually measured experimentally us- 
ing the ESS ansatz O], i.e. via the relative scaling 
(\A T v\i) ~ (|A T v| 2 )^Vwhere £(q) = C(?)/C(2). In Kol- 
mogorov K41 phenomenology 0, the Lagrangian second 
order structure function is assumed to scale linearly with 
the time increment r so that £(2) is considered as being 
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FIG. 2: PDF n T (Au) calculated for time lags r = 
1.3,2.7,5.4,11.2,22.4,44,89.3 and 174% of T L . The curves 
are computed for variations normalized to the standard de- 
viation and displayed with vertical shift for clarity. The un- 
derlying crosses in the outer curve correspond to the MRW 
model — see text — with A 2 = 0.115. 

equal to one, eventhough it is quite hard to establish it 
experimentally 0. Given the statistics available in our 
experiments, one is limited to moments up to order 6. 
We observe that £(g) is well represented by a quadratic 
law 0: 

Z(q) = (l/2 + \l)q-\lq 2 /2, Xl = 0.115 ± 0.01 . (1) 

Note that the Lagrangian value of the intermittency pa- 
rameter (\\ = 0.115) is larger than the value (A|; = 
0.025), measured from hot-wire anemometry or in direct 
numerical simulations |L2j. This is expected [jl3|, and 
mainly comes from the facts that time increments are 
measured here and that the reference structure function 
is the second order one. 

Intermittency is thus observed and quantified in both 
Lagrangian and Eulcrian frameworks. In contrast to 
traditionnal Eulerian studies where intermittency is de- 
scribed in terms of multiplicative processes, we look here 
for a dynamical origin. We consider the statistics of the 
fluid particles fluctuating velocity in analogy with a ran- 
dom walk. We write a velocity increment over a time lag 
r as the sum of contributions over small times t\ : 

r/n 

A T v(t) = v(t + r) - v(t) = ^2 At!«(< + ran) • (2) 

n=l 

If the incremental "steps" of duration n were inde- 
pendent (and identically distributed), the PDF H T (Av) 
would readily be obtained as a convolution of the ele- 
mentary distribution at scale T\ , n Tl (Av) — plus an 
eventual convolution kernel to account for stationarity 
at large scales. Such a regular convolution process corre- 
sponds to the Kolmogorov K41 picture of turbulence [0; 
the particle velocity fluctuations are Brownian and the 
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scaling is monofractal. A first important result of our 
analysis is to show that the elementary steps are not inde- 
pendent. The auto-correlation of the signed increments, 
A T1 v(t), decays very rapidly (cf. Fig. 3): the correlation 
coefficient drops under 0.05 for time separations larger 
than 2ti. However, if one considers the amplitude of the 
"steps" (\A Tl v(t)\), one finds that the auto-correlation 
decays very slowly and only vanishes at the largest time 
scales of the turbulent motion. Recast in terms of the 
random walk, our results show that the amplitudes of the 
"steps" are long-range correlated in time although their 
directions are not. As this point is fundamental for our 
approach, we have verified it using a Lagrangian tracking 
algorithm in a Direct Numerical Simulation (DNS) of the 
Navier-Stokes equations, using a pseudo-spectral solver, 
at Ri — 75, and for the same ratio Ti/Tj, — see inset of 
Fig. 3. The results are in remarkable agreement with our 
measurements. All increments are correlated for At < t\, 
the time over which they are computed. Above t±, the 
correlation of the signed increments rapidly drops while 
the correlation coefficient of their absolute values decays 
very slowly, to vanish only for At > 3 Tl ■ This behavior 
is observed for t\ chosen from the smallest resolved time 
scales to inertial range values. These observations also 
persist in the limit of very small time increments, and 
thus presumably for the acceleration of the fluid particle 
and thus for the forces acting on it. 
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FIG. 3: Variation of the normalized correlation coeffi- 
cient x(f,9)(At) = ((/(* + At) - (f))(g(t) - (g)))/a f a g . 
Two velocity components are considered: the squares and 
diamonds mark the x{A T1 v x , A T1 v x ) and x(A T iVy, A T1 v y ) 
auto-correlation functions while the circles mark the cross- 
correlation x(A T1 Vx, A T1 v y ); curves with filled symbols are 
computed using the absolute value of the increments while 
the curves with open symbols are computed using the full 
signed increments. The main curve corresponds to the ex- 
periment at Ri = 740, with n = 0.03 Tl- The inset shows 
similar results for a Direct Numerical Simulation at Re = 75. 

Theoretically, one would like to understand this behav- 
ior from the hydrodynamic forces in the Navier-Stokes 
equations. Such a direct analytical treatment is out 
of reach at present. We propose, as a natural first 



step, to study a surrogate dynamical equation of the 
Langevin type. In this procedure, one considers a one- 
dimensional variable, W(t), representing the particle ve- 
locity, driven by a stochastic force. If this force is chosen 
as a white noise then W(t) has the dynamics of Brow- 
nian motion: its statistics is monofractal with a simi- 
larity exponent equal to 1/2 — the increments scale as 
(\W(t + t) — W(t)\ p ) t ~ r p / 2 , corresponding to the non- 
intermittent Kolmogorov 1941 picture. In order to ac- 
count for intermittency, one needs to ascribe other prop- 
erties to the stochastic force. Guided by our experimen- 
tal results, we build a stochastic force having a random 
direction and a long-range correlation in its magnitude. 
Specifically, its direction is modeled by a Gaussian vari- 
able G(£), chosen white in time, with zero mean and unit 
variance. The amplitude of the force, A(t), being a posi- 
tive variable, is written A(t) = exp[uj(t)] where the mag- 
nitude ui(t) is a stochastic process that satisfies: 

(uj(t)uj(t + At)) t = -A 2 ln(At/T L ) for At < T L , (3) 

and otherwise — A 2 being an adjustable parame- 
ter. When discretized, this dynamics corresponds to a 
one-dimensional Multifractal Random Walk (MRW) Q . 
Analytical calculations show that the resulting dynam- 
ical variable W(t) has multi-scaling properties. The 
moments have scaling laws, (|A T IT| 9 ) ~ r^ q \ with 
C(q) = (1/2 + A 2 )<? - AV/2, so that A 2 in equation (g) 
is the intermittency parameter of the model juj. It is a 
fundamental point that the same parameter A 2 governs 
both the evolution of the PDFs of the increments (one- 
time statistics) and the time correlation of the process 
(two-time statistics). 

We show that this model captures the essential features 
of the Lagrangian data. First, in order to test the rele- 
vance of equation |[ we have computed, from the experi- 
mental and numerical data, the auto-correlation function 
of the logarithm of the amplitude of infinitesimal La- 
grangian velocity increments: x(log | A Tl v x | , log | A Tl v x \). 
Figure 4(a) confirms that the logarithmic decrease build 
in the MRW model (equation g) is observed both in 
the experimental and numerical data; it yields the es- 
timate A 2 = 0.115 ± 0.01. Second, we check the rele- 
vance of the model for the description of the one-time 
statistics of the Lagrangian increments A T v. We note — 
Fig. 2, upper curve — that the choice A 2 = 0.115 yields 
a PDF for the stochastic force that is in remarquable 
agreement with experimental measurements of fluid par- 
ticle accelerations 0j. The agreement at larger time 
scales is evidenced on the behavior of the first two cu- 
mulants. Cumulants are computed with more reliabil- 
ity that the moments and are related to them through 

(\A T vn = (exp(qln|A rV |)} =exp(E„C„(r)qV?!)- In 
the MRW model, one can analytically derive |]l4j| : 

C 1 (r) = (l + A 2 )ln(r), C 2 (t) = —A 2 ln(r), (4) 
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dominate the motion of a fluid particle. It can be un- 
derstood if, along its trajectory, the particle encounters 
very intense small-scale structures (vortices and stagna- 
tion points) over a more quiet background. Intermittency 
is then due to the nature and distribution of these small 
scale structures. Indeed, the analogy with a random walk 
suggests that the statistics at all scales can be recov- 
ered if one ascribes two properties to the small scales: 1) 
the probability density function of fluid particle acceler- 
ations and 2) the functional form of their time correla- 
tions. In the Lagrangian framework, these features are 
directly linked to the Navier-Stokes equations that gov- 
ern the elementary changes in the velocity (momentum) 
of the fluid particles. It thus gives a possibility to derive 
intermittency from the constitutive physical equations. 
Although this may be quite a theoretical challenge, di- 
rect numerical simulations look promising as they allow 
the study of the flow dynamical fields (pressure, velocity 
gradient tensor, etc.) along the trajectory of individual 
fluid particles. 
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FIG. 4: Experimental (open squares) and numerical data 
(open triangles) compared to the predictions of the MRW 
model (filled circles). (a) Correlation in time of the 
magnitude of one component of the velocity increments, 
x(log |A T1 «.c|,log IAtjUjuQ, computed for a time lag ri = 
0.03Tl. (b) & (c): first and second order cumulants versus 
time scale r. 

all higher order cumulants being null. Ci(r) and C2(t) 
computed from the experimental and numerical data are 
shown in Figures 4(b) and 4(c) and compared to MRW 
model predictions when the intermittency parameter is 
set to the value A 2 = 0.115 that is derived from the cor- 
relations in the dynamics. One observes that in each 
case the agreement is excellent; the slope of the variation 
9Ci } 2(t) /<91nr in the inertial range is correctly given by 
equation (|J). The same intermittency parameter thus 
governs the anomalous scaling of the Lagrangian velocity 
increments and their long-time dynamical correlations. 

We therefore believe that long-time correlations in the 
Lagrangian dynamics are a key feature for the under- 
standing of intermittency, which leads to a new dynam- 
ical picture of turbulence. Long-time correlations and 
the occurrence of very large fluctuations at small-scales 
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